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Abstract Our recent studies of pulsar population statistics suggest that 



improvements of radio and gamma-ray beam geometry and luminosity mod- 
els require further refinement. The goal of this project is to constrain the 
viewing geometry for some radio pulsars, especially three-peaked pulse pro- 
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files, in order to limit the uncertainty of the magnetic inclination and impact 
angles. We perform fits of the pulse profile and position angle sweep of ra- 
dio pulsars for the available frequencies. We assume a single core and conal 
beams described by Gaussians. We incorporate three different size cones 
with frequency dependence from the work of Mitra & Deshpande (1999). 
We obtain separate spectral indices for the core and cone beams and explore 
the trends of the ratio of core to cone peak fluxes. This ratio is observed to 
have some dependence with period. However, we cannot establish the sug- 
gested functional form of this ratio as indicated by the work of Arzoumanian, 
Chernoff & Cordes (120021) . 
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1 INTRODUCTION 

A recent population synthesis study (Gonthier et al. I2004[) that incorporates the radio 
beam geometry and luminosity models of Arzoumanian, Chernoff & Cordes (2002) (here- 
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after ACC) and the geometry and luminosity of the polar cap slot gap model of Muslimov 
& Harding ( 2003) simulates the number of radio-loud and radio-quiet (with a radio flux 
below the survey threshold) 7-ray pulsars detected by the instruments EGRET, AGILE 
and GLAST. The simulated correlations between the radio and bright 7-ray beam pro- 
files suggest that 7-ray profiles with two peaks have a radio profile that is core dominated 
with the core peak appearing in between the two 7-ray peaks, which is contrary to the 
features of most of the EGRET detected 7-ray pulsars that exhibit a single radio peak 
leading in time the two 7-ray peaks (Thompson I200lj) . Core-dominated radio emission 
is implied in the ACC model, especially for young 7-ray pulsars, due to the assumption 
that the radio core- to-cone peak flux has a 1/P period dependence. These findings sug- 
gest that a better radio beam model is needed to account for the observed correlations 
of radio and 7-ray profiles. Motivated by these results, we seek to develop an alternate 
radio beam geometry and luminosity model by focusing on well defined three-peak ra- 
dio pulse profiles with adequate polarization data. We fit the polarization position angle 
within the rotating vector model (RVM) of Radhakrishnan & Cooke (1969) as well as 
the pulse profiles. While we realize that other possible interpretations of the structure of 
radio beams are possible, such as the core being a conal structure (Kijak & Gil 2003]) or 
the emitting surface being patchy (Han & Manchester 2001]) , we assume that pulsars are 
standard candles, there is a single core beam and a single cone beam characterized by 
an inner, middle or outer width following Mitra & Deshpandc (1999), and the rotating 
vector model is applicable. 

2 MODEL 

The angular distribution of the core beam is assumed to be a Gaussian centered along 
the magnetic axis with a characteristic width represented by 



where P is the period in seconds and the factor in front of the exponent normalizes the 
Gaussian. For a given viewing geometry, defined by the inclination angle, a, from the 
magnetic axis and the viewing angle, C, relative to the rotational axis, then the polar 
angle, 9, relative to the magnetic axis is related to the phase angle, 0, which is the 
azimuthal angle relative to the rotational axis, through the expression 



The 4>o° le is the phase offset associated with the core beam of the pulsar profile. The 
characteristic width, p C orc, was adopted from ACC model. 
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cos 9 = sin a sin £ cos(4> — 4> ) + cos a cos £■ 



(2) 
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The angular distribution of the emitting surface for the conal beam is also parame- 
terized with a Gaussian by the expressions 

W<U) = 9 3/2 m T7~\ l^' ati ° - !) f 1 - sin (? sin ^)l + 2 | e-V-WW, where 

4 °. 8 (^1 + ^P- 1/2 , 



= ^ r p CO no, and 

^w/^cone 



4Vln2 

where £ is the azimuthal angle in the magnetic frame and v is the observing frequency 
in MHz. The fitting parameter /3 ra tio is introduced to describe an asymmetric conal 
distribution in order to have the ability to fit conal peaks in the pulse profile that exhibit 
different maxima. The two fitting parameters, S r and S w , allow for various cones with 
(5 r =0.8, 1.0, and 1.3 for the inner, middle, and outer cones of Mitra & Deshpande (|1999p . 
While we do not believe that the conal beam actually has this particular geometry as 
opposed to say a patchy distribution, this function merely allows us to describe the heights 
of the two conal peaks in the pulse profile with a function that is easily integrated and, 
therefore, normalized. 

We fit the position angle, ip, of the polarization vector with the RVM using the 
relationship 

tan^ - = . - (4) 

sin (, cos a — cos C sin a cos(<p — <Po ) 

The phase offset angle, (f>g A and the position angle offset, tp , can be obtained for each of 
the available frequencies, then subtracted from the originals to obtain a global position 
angle sweep for all frequencies. The inclination angle a and the impact angle (3 are related 
to each other by the maximum rate of position angle (PA) sweep V4iax> 

(dip\ sin a 

'/'max = "77 = 77777 • (5) 

V#/max sm P 

Without high quality polarization data, it is not possible to uniquely determine the 
viewing geometry, a and £, from the position angle alone. We fit the global position 
angle to obtain the maximum rate of change of the position angle, V'max- While the Stokes 
parameters I, Q, U and V have standard deviations that obey Gaussian statistics, values 
derived form these parameters do not. We follow the work of Everett & Weisberg (2001) 
in assigning errors to the position angle, ip. Presented in Figure 1 is the global fit of the 
PA for pulsar, B2045-16. 

Assuming the maximum rate of change of the position angle V'maxj we u ^ t ne profiles 
constraining the inclination angle a, and obtain the impact angle, [3. Errors in the profile 
are determined by the standard deviation of the off beam noise. The fitting parameters 
associated with a component of the profile are the overall amplitude, Aj, and the phase 
offset, where z=core or cone. The flux in the profile is given by the expression 

s i (6,v)=A i {u)f i {6,v), (6) 



(3) 
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Fig. 1 The global position angle of the pulsar, B2045-16, that combines the 
data obtained of frequencies of 430, 606, 925, 1408 and 1624 MHz. The smooth 
curve is the RVM fit. 

where the /j term is defined above in Equations 1 and 3. Since the angular distributions 
are normalized over the emitting surface, the coefficients, Ai, represent the total angle 
integrated core and cone flux at the particular frequency, v. However, the width, 5 W , 
and radius, <5 r , of the conal annulus are allowed to vary to account for the various cones 
suggested by the observations. In addition, the parameter, /3 ra tio allows for the two peaks 
of the conal beam in the pulse profile to have different contributions. Assuming a power 
law with a low frequency cutoff of 50 MHz, the frequency-differential flux has the form 



where i = core or cone, is the spectral index, d, is the assumed distance to the pulsars 
and Li is the angle and frequency integrated luminosity. Having the amplitudes as a 
function of the available frequencies, the spectrum can be fit with this function to obtain 
the luminosity and spectral index of the core and conal components, separately. 

Since the maximum rate of change of the PA, tp'mux' * s the slope of the sweep at 
the inflection point of the curve, the minimum \ 2 of the RVM fit to the PA is weakly 
dependent on the inclination angle, a, as seen in Figure 2(a), where the \ 2 °f the fit is 
plotted as a function of the maximum rate of change in the PA, ip' max , and a. On the 
other hand, the \ 2 obtained from fitting the profile is weakly dependent on the maximum 
rate of change of the PA, ip' max , as seen m Figure 2(b). Therefore, fitting both the PA 
and the profile gives a better constraint of the viewing geometry. 

The largest uncertainty of the procedure is introduced by the selection of the radius 
of the annulus parameter, <5 r , which defines the inner, middle or outer cone. Values of 
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Fig. 2 The reduced x 2 obtained from fitting the PA (a) and the profile (b) of 
the pulsar B2045-f6 as a function of the maximum rate of change of the PA, 
Vmax' ano - t ne inclination angle, a, at a frequency of 430 MHz. 

S T = 0.8, 1.0 or 1.3 result in different inclination angles. For example in the case of the 
pulsar, B2045-16, once the maximum rate of change of the PA is established and allowing 
all parameters to vary while fitting the profile, the minimum \ 2 — 294 is obtained for 
S t = 1.00, <5 W = 0.65 and a = 27°. 9, suggesting the preference for the middle cone. 
However, as can be seen in Table 1, there is a large variation in the values of a with the 
conal radius parameter, 6 T . Based on the average minimum \ 2 of 265, the middle cone is 



Table 1 Averages and standard deviations of the parameters obtained from 
fitting the profiles of the pulsar, B2045-16, for the inner, middle and outer cones 
at the frequencies of 408, 606, 925, 1408 and 1642 MHz. 



Cone 


5r 


5 W 


/3ratio 


a 





x 2 


inner 


0.8 


0.50 ±0.03 


0.76 ±0.13 


21.7 ±0.6° 


-0.67 ±0.02° 


356 ± 410 


middle 


1.0 


0.65 ±0.04 


0.69 ±0.10 


27.8 ± 0.7° 


-0.84 ±0.03° 


265 ± 293 


outer 


1.3 


0.93 ±0.08 


0.65 ±0.09 


37.6 ± 0.9° 


-1.10 ±0.04° 


544 ± 706 



selected, and we show in Figure 3 the resulting profiles and fits of the pulsar, B2045-16, 
for the indicated frequencies. 

In Figure 4, we plot the core and cone coefficients, ^4 C orc and A corLei representing the 
angle integrated flux of Equation 6 from the fits of the profiles of pulsar, B2045-16, as a 
function of the available frequencies indicated. The smooth curves are fits of the power 
law given in Equation 7, resulting in the indicated total angle and frequency integrated 
fluxes and spectral indices for the core and cone beams. Using the distances (with «20% 
uncertainty) of the pulsars from the new distance model of Cordes & Lazio p002[) . we 
can estimate the core and cone luminosities. 
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Fig. 3 The pulse profiles (curves with dots) and the fits (smooth curves) of the 
pulsar B2045-16 at the indicated frequencies. 

3 RESULTS 

We selected a group of seventeen pulsars ranging in period from 0.06 s (B1913+16) to 
3.7 s (B0525+21) whose pulse profiles exhibit a fairly well defined structure consisting 
of three peaks and have adequate polarization data. We obtained the data from the 
European Pulsar Database (EPN) with all profiles and polarization data taken in the 
survey by Gould & Lyne (|1998p . From the fitting of the position angle sweep and pulse 
profiles, we obtain the inclination angles, a, the impact angle, /3, the spectral indices, 
fluxes and luminosities of the core and cone beams separately. In Table 2, we compare 
the extracted inclination a, impact angles (3 and the maximum rate of change of the 
PA, VCax> witn tne works of Lyne & Manchester (|1988)) and Rankin (|1993j) . With this 
information, we are able to estimate the ratio of the core-to-conc peak fluxes, which is 
independent of the pulsar distance and associated uncertainties. In Figure 5, we present 
this ratio as a function of the period for three different indicated frequencies. The pre- 
dictions of the ACC model that assume the ratio of core-to-cone peak flux has a 1/P 
dependence are indicated by dotted lines in each of the panels. The period dependence of 
the ratios indicates a break around a period of 0.7 seconds. We parameterize the period 
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Fig. 4 The core and cone angle integrated fluxes obtained from the fits to the 
profiles of the pulsar B2045-16 at the frequencies of 408, 606, 925, 1408, and 
1642 MHz. Assuming the power law in Equation 7, the fluxes are fit to obtain 
the indicated angle and frequency integrated fluxes and spectral indices. 



and frequency dependence of the ratios below and above 0.7 s. We find that the ratio of 
the core-to-cone increases as P 15 v~ 1 up to 0.7 seconds and then decreases as P~ z v~ x . 
The results are quite different than those of the ACC model, particularly for the behavior 
of the ratio core-to-cone peak flux for short period pulsars. Our results suggest that they 
are not as core-dominated as implied by an extrapolation of the ACC model. 

With the extracted fit parameters, we are able to estimate the luminosity of the core 
and cone beams separately, as shown in Figure 6. The bottom part of the figures displays 
the luminosity of the cone beam (solid dots) with a parameterized function suggested 
by the fits (crossed circles). The functional dependence of period and period derivative 
is indicated in the figure for the parameterization. Likewise the luminosity of the core is 
displayed in the middle portion of the figure. The total luminosity is shown at the top of 
the figure along with the predictions of the ACC model. 

4 CONCLUSIONS 

We focus on a select group of pulsars that exhibit three peaks in their pulse profiles with 
adequate polarization data. We find that there are a very limited number of these pulsars 
on the EPN database with good polarization data. All of polarization data in this study 
comes from the survey of Gould & Lyne (|1998|) . We fit both the position angle as well as 
the pulse profile intensities with a simple model that assumes a single core beam and a 
single cone beam. Assuming a standard period dependence of the characteristic widths 
(Equations 1 and 3) and the presence of three conal beams as described by Mitra & 
Deshpande (|1999[) . we find that the ratio of the core-to-cone peak flux for short period 
pulsars does not follow the 1/P trend suggested by the study of Arzoumanian, Chcrnoff 
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Table 2 Inclination angles a, impact angles /3, and maximum rate of change of 
PA ip' max obtained from fits are compared to the studies of Lyne & Manchester 
(T19881) and Rankin (|T993|l . In the fits a < 90° and (3 and ip' max have the same 
signs. 



Pulsar 




This Work 


±jy lie 


& Manchester 




Rankin 


Name 


a 


R 

H 


ib' 


Q 




ib' 


a 


3 


ib' 

rmax 


B1913+16 


44 


-0.8 


-52 








46 


-0.8 


-51 


B 1804-08 


51 


3.5 


13 


47 


2 


21 


63 


5.1 


10 


B1702-19 


45 


5.7 


7 


90 


4.1 


14 


85 


4.1 


14 


B2048-72 


35 


-2.3 


-14 


29 


-1.4 


-20 








B1839+09 


90 


2.6 


22 


90 


2.9 


20 


83 


1.4 


42 


B0727-18 


25 


-4 


-7 


28 


-6.8 


-4 








B2003-08 


19 


-3.7 


-5 


13 


-3.1 


-4 


13 


-3.3 


-4 


B1845-01 


42 


-3.2 


-12 














B1508+55 


35 


2.9 


11 


80 


2 


28 


45 


2.7 


15 


B1821+05 


27 


2.9 


10 


28 


1.5 


18 


32 


1.7 


18 


B2111+46 


10 


1.1 


9.1 


8.6 


1.3 


6.7 


9 


1.4 


6.7 


B1039-19 


24 


1.3 


18 


34 


1.8 


18 


31 


1.7 


18 


B2045-16 


27 


0.8 


31 


37 


1.1 


30 


36 


1.1 


30 


B0525+21 


14 


-0.4 


-35 


23 


-0.7 


-31 


21 


-0.6 


-36 



& Cordes (2002). There appears to be a break in the period dependence of the ratio near 
a period of 0.7 seconds, with the ratio decreasing for short period pulsars. Such a result 
alters significantly the interpretation of radio profiles of short period pulsars and requires 
further study. 

The interpretation of radio profiles in young high-energy pulsars is not well estab- 
lished. While Arzoumanian, Chernoff & Cordes (2002) proposed both a beam geometry 
and a luminosity model of the core and cone radio beams, a considerable controversy 
exists in the interpretation of single peaked radio pulse profiles, as to whether the origin 
is from core emission or partial cones (Manchester I2001|) . The study presented in this 
paper and others like it may help in understanding the geometry of radio emission from 
pulsars, particularly from young pulsars. Our study indicates that the ACC model re- 
quires some revision in regards to the ratio of the core-to-cone emission for short period 
pulsars. 
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Fig. 5 The logarithm of the ratio of the core-to-cone peak flux is plotted as 
a function of the logarithm of the pulsar period for the indicated frequencies. 
The predictions of the ACC model are indicated by dotted lines. The dashed 
and solid lines represents the parameterization of the ratios below and above 
0.7 seconds, respectively (see text). 
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Fig. 6 The logarithm of the total (top), core (middle) and cone (bottom) lu- 
minosities obtained from the fits (solid dots), parameterization (circles with 
crosses) and predictions of ACC model (top). 
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